GAMOW-TELLER STRENGTH DISTRIBUTIONS 

for /3/3-DECAYING NUCLEI WITHIN 

CONTINUUM-QRPA 

■^ ; S.Yu. Igashov\ V.A. Rodin^ M.H. Urin^'^ A. Faessler^ 

^ ; ^ Moscow Engineering Physics Institute (State University), Russia 

CN i ^Institute for Theoretical Physics, University of Tuebingen, Germany 

^ '. ^Kernfysisch Versneller Instituut, University of Groningen, The Netherlands 

:z; 

OS 



a 



o 



February 2, 2008 



Abstract 

Q ■ A version of the pn-continuum-QRPA is outlined and applied to describe the 

2 . Gamow- Teller strength distributions for /3/3-decaying open-shell nuclei. The calcu- 

lation results obtained for the pairs of nuclei ^^^Cd-Sn and ^^'^Te-Xe are compared 
with available experimental data. 

PACS numbers: 25.40.Kv, 23.40.-s, 24.30. Gd, 27.50.+e 



in 

^ ■ 1 Introduction 



Description of weak interaction in nuclei is often a challenge for models of nuclear struc- 
ture. Numerous calculations of the nuclear /3/3-decay amphtudes well illustrate this state- 
ment (see, e.g. Refs. [1] and references therein). Uncertainties in theoretical calculations 
of the Gamow- Teller (GT) 2z//5/3-decay amplitude M^^ have stimulated experimental 

rS I studies of the GT^^^-strengths of the 1"*" states virtually excited in the decay process (see, 

_c^; e.g. Refs. [21 [3]). 

As a double charge-exchange process, 2z/(0i/)/3/3-decay is enhanced by nucleon pairing 
which is due to the singlet part of the particle-particle (p-p) interaction. The discrete 
quasiboson version of the quasiparticle RPA (pn-dQRPA) which accounts for the nucleon 
pairing is usually applied to calculate the /3/?-decay amplitudes in open-shell nuclei p. 
In spite of differences in model parameterizations of the nuclear mean field and residual 
interaction in the particle-hole (p-h) and p-p channels, all pn-dQRPA calculations reveal 
marked sensitivity of the amplitude M|.^ to the ratio gpp of the triplet to singlet strength 
of the p-p interaction. Physical reasons for such a general feature of all calculations were 
analyzed in Ref. [1], where they were attributed to violation of the spin-isospin SU(4) 
symmetry in nuclei. An identity transformation of the amplitude into sum of two terms 
was used in Ref. [5|. One term, which is due to the p-p interaction only, depends linearly 
on gpp and vanishes at gpp = 1 when the SU(4)-symmetry is restored in the p-p sector 
of a model Hamiltonian. The second term is a smoother function of gpp at gpp ~ 1, but 



exhibits a quadratic dependence on the strength of the mean-field spin-orbit term, which 
is the main source of violation of the spin-isospin SU(4)-symmetry in nuclei. 

Understanding of general properties of the amplitude M^.^ helps to improve reliability 
of evaluation of /3/5-decay amplitudes. For a quantitative analysis, we use here an isospin- 
selfconsistent pn-continuum-QRPA (pn-cQRPA) approach of Ref. [5], where this approach 
was applied to evaluate the GT*^^) strength distributions in single-open-shell nuclei. In 
the reference the full basis of the single-particle (s-p) states was used in the p-h channel 
along with the Landau-Migdal forces, while the nucleon pairing was described within 
the simplest version of the BCS-model based on discrete basis of s-p states. A rather 
old version of the phenomenological isoscalar nuclear mean field (including the spin-orbit 
term) was used in Refs. [5], H], as well. 

The first application of the pn-cQRPA approach of Ref. [5] to description of the 
/3/3-decay observables in several nuclei has been given recently in Ref. [7j. Realistic (zero- 
range) forces have been used in the p-p channel to describe the nucleon pairing within 
the BCS model realized on a rather large discrete-fquasidiscrete s-p basis. 

The pn-cQRPA approach of Refs. [H [7| is further extended here by using a modern ver- 
sion of the phenomenological isoscalar mean field (including the spin-orbit term) deduced 
in Ref. |6] from the isospin-selfconsistent analysis of experimental single-quasiparticle 
spectra in double-closed-shell nuclei. In the present contribution we give a brief overview 
of the approach and its applications to description of different GT strength functions for 
the pairs of nuclei ^^^Cd-Sn and ^^°Te-Xe. 

2 Coordinate representation of the pn-cQRPA equa- 
tions and GT strength functions 

In formulation of a version of the pn-cQRPA we follow Ref. [5], where the pn-dQRPA 
equations originally written in terms of the forward X^~^ and backward amplitudes YJ~^ 
are transformed in equivalent equations for the 4-component radial transition density 
Pj (s, r). The latter is defined in terms of the X, Y amplitudes by Eqs. (39), (40) of Ref. 
[S] and related to the GT^"-* excitations having the wave functions |l+/i,s) and energies 
LUg. The expression for the transition density p\ {s, r) related to GT'^^^ excitations follows 
from that for p~ [s^r] by the substitution p ^^ n. Hereafter, we often use the notations 
of Ref. [S] (except for the energies) and refer to some equations from this reference. 
Limiting ourselves in this contribution to description of the GT transitions only, we omit 
the quantum numbers indices J = S = 1, L = for spin-monopole excitations. The 
spin-angular variables in all expressions are separated out as well. 

The pn-dQRPA solutions cUs are related to the excitation energies EJ^J measured from 
the ground-state energy Eq{Z ±1,N ^ 1) of the corresponding daughter nuclei as: 

^. ± (/ip - /in) = ^f ) = i^g) + Q^^ (1) 

Here, yUp(„) is the chemical potential for the proton (neutron) subsystem found from the 
known BCS equations, Q^^ = £b{Z, N) — £b{Z ± 1, A^ =i= 1) are the total binding-energy 
differences, cj^^^ are the excitation energies measured from Eq(Z, N) —Y.a i^^a = ~£biZ, N) 
{iria is the nucleon mass). The energies cj^^ are usually described by a model Hamiltonian. 



The system of equations for pf (s, r) (Eq. (41) of Ref. [5]) contains the exphcit expres- 
sion for the 4x4 matrix of the free two-quasiparticle propagator A^-^ (r, r', Ug) (Eq.(43) of 
Ref. [5j; Ai^^ = Al/. (p •«-»• n)). These propagators are the main quantities in description 
of charge-exchange excitations within the pn-QRPA. In particular, in terms of Aik one 
can formulate a Bethe-Salpeter-type equation for the effective propagator Aik{r,r',Us) 
[?]. The spectral expansion of Aik in terms of pi{s, r) allows one to express the pn-QRPA 
strength functions in terms of the effective propagator, or, equivalently, in terms of the 
4-component effective fields [5]. Some relevant formulas are shown below. 

The GT^^^ strength functions, corresponding to the external fields (probing operators) 
V;^^^ = Ea V^^^a), V;(^^ = (T^,T^^\ are defined as follows: 

S^^\uj) = j:\{l\s\\V^^^\\0)\'6{uj-ujg) (2) 



with GT strengths Bi^^^GT) = |(1+, s||yW||0)p. The strength function S^-\uj) can be 
expressed in terms of the corresponding effective field V^L (r, a;), which is different from 
the external one V^ (r) = 5ii due to the residual interaction [5]: 

S^'\uj) = --ImY, f A[-;\r,r',uj)v}^;^\r',u;)drdr', (3) 

The residual interaction here is supposed to be of zero-range type with intensities F^ . 
Fi = F2 = 2G', -F3 = F4 = Gi. For the O"*" p-h and p-p channels the corresponding 
strengths are: F{ = Fg = 2F' and F^ = F^ = Gq, respectively. Dimensionless values 
g' = G'/C, f = F'/G, {G = 300 MeV ■ fm^) are the well-known Landau-Migdal p-h 
strength parameters. The same parameterization we use for the p-p interaction strengths: 
91 = Gi/G, go = Gq/G. For calculation of S^~^\uj) one can use Eqs. (3), (4) with 
substitution p ^->- n [5J. An alternative way is based on the symmetry properties of Aik'. 
All = ^22 • As a result, we get the expression for S^'^^cu) in terms of Ai^ and V^ . 
This expression is obtained from Eqs. (3), (4) with the substitution 1^2. 

The nuclear GT*^"^ amplitude for 2z//3/3-decay into the ground state |0') of the product 
nucleus (A^ — 2, Z -|- 2) is given by the expression: 

Mg^-E^"''""''"''''''''-'"^"""'. (5) 

s ^s 

where ujs = E^ - |(Fo + ^o') = E^,s + ^{Qb^^ + Qb^^')- To calculate M^^ within the 
pn-QRPA, the vacua |0) and |0') should be identified. As a result of such identification, 
one has COs = ^{^i~^ + ^i'^^ ) ~ ^s, in accordance with Eq. (1). 

The amplitude (5) can be expressed in terms of a "non-diagonal" GT*^^) strength 
function S^ Huj): 



M^^y = / uj-^S^—\u)du, (6) 

where S^ \uj) is defined as follows: 

S^--\u;) = Y.{0'\\V^~^\l+,s){l^,s\\V^~^\\0)6{uj - Us). (7) 



The corresponding pn-QRPA expression for S^ ^ is: 

S^—\uj) = --ImJ2 f A^\r,r',uj)V.[;^\r',ij)drdr'. (8) 

i 

An alternative expression for M^^ is obtained in terms of the "non-diagonal" static 
polarizibility [7]: 



M'Jt 



-\Y.J Ai\^y^^ = 0)V^iV(^',^ = ^)drdT'. (9) 



Decomposition of the amplitude (5) into two terms [1] 

M^^ = {Ml^r^y + u}^^T.j^EWSR^—\ (10) 

i?iy5i?(-) = ^cD,(0'||i>(-)||l+,s)(l+,s||v(-)||0), (11) 

s 

where ujgtr is the energy of GT*^^) giant resonance (GTR), allows us to clarify the sensitive 
dependence of 2i//3/3-decay amplitude as a function of Qpp (for details, see Ref. [1]). The 
"non-diagonal" energy-weighted sum rule EWSR^ ^ is straightforwardly expressed in 
terms of the strength function S*^ -* of Eq. (6): 

EWSR^--^ = f ujS^—\u)duj, (12) 

again supposing the QRPA vacuum |0') is identified with that of |0). 

3 Calculation of strength function within the pn-cQRPA 

Starting from the coordinate representation of the pn-dQRPA equations outlined above, 
we are able to take exactly into account the s-p continuum in the p-h channel and, there- 
fore, to formulate a version of the pn-cQRPA. The pairing problem is solved on a rather 
large basis of bound-|-quasibound proton and neutron s-p states within the present ver- 
sion of the model. To take the s-p continuum into account, the following transformations 
of the expression for Aii:{r,r',uj) [5] are done: (i) the Bogolyubov coefficients vx, u\ 
and the quasi-particle energies Ex are approximated by their non-pairing values v\ = 0, 
ux = 1, and Ex = ex — fi for those s-p states (A), which lie far above the chemical po- 
tential (i.e. Ex — H^ Aa), (ii) the Green function of the s-p radial Schrodinger equation 
5'(A)(r,r',e) = I]e^(£ — £a + ^0)^^xa(^)Xa('^')) "which is calculated via the regular and irreg- 
ular solutions of this equation, is used to perform explicitly the sum over the s-p states in 
the continuum. As a result, the properly transformed free two-quasiparticle propagator 
A is obtained, upon which a corresponding version of the pn-cQRPA is based. 

The solution of the pairing problem is simplified by using the "diagonal" approxi- 
mation for the p-p interaction for the 0^ neutral channel. In this approximation the 
nucleon-pair operators are assumed to be formed only from the pair of nucleons occupy- 
ing the same s-p level A. The nucleon pairing is described with the use of the Bogolybov 
transformation with the gap parameter Aa dependent on A. The same number Nb+qt of 
bound+quasibound states forming the basis of the BCS problem is used for both the neu- 
tron and proton subsystems. These numbers are shown in Table 1 for nuclei in question. 



In evaluation of total binding energies within the model (that is necessary to evaluate the 
pairing energies Epair) the blocking effect for odd nuclei is taken into account. In descrip- 
tion of the nucleon pairing, different values of the p-p interaction strength parameters 
QQ^n and Qq^p for the neutron and proton subsystems are used. These values are found 
from comparison of the calculated and experimental pairing energies for nuclei under 
consideration (Table 1). 



Table 1: The phenomenological mean field parameters (f/o, Uso and a), singlet and triplet 
p-h and p-p interaction strengths (/', g' , go, g^p) used in calculations. 



Pair of nuclei 


Uq, MeV 


t/so, MeV-fm2 


a, fm 


/' 


90, n 


9o,p 


Nb+qb 


9' 


9pp 


ii6cd-ii6Sn 


51.62 


34.08 


0.618 


1.06 


0.388 


0.333 


22 


0.77 


1.0 


i30xe-i30xe 


51.74 


34.025 


0.628 


1.09 


0.356 


0.364 


22 


0.88 


0.99 



The mean field consists of the phenomenological isoscalar part (including the spin- 
orbit term) along with the isovector and Coulomb part (Eq. (1) of Ref. [5]). The 
parameterization of the Woods-Saxon-type isoscalar part contains two strength (f/o, Uso) 
and two geometrical (rg, a) parameters [6J. The mean field isovector part (the symmetry 
potential) is calculated in an isospin-selfconsistent way (Eqs. (7), (35) of Ref. [5]) via the 
neutron-excess density and Landau-Migdal strength parameter /'. The mean Coulomb 
field is also calculated selfconsistently via the proton density. All densities are calculated 
with taking into account the nucleon pairing. Five above-listed model parameters found 
in Ref. [6] for a number of double-closed-shell nuclei are properly interpolated for nuclei 
under consideration (see Table 1; tq = 1.27 fm is taken for all nuclei). 

The values of the Landau-Migdal strength g' listed in Table 1 are obtained by fitting 
the experimental GTR energy in calculations of the GT*^^) strength function. The p- 
p interaction strength gi (or its relative value gpp = 2gi/{go^n + 9o,p)) is considered as 
a free parameter. It can be adjusted to reproduce the experimental M^^ value (the 
corresponding values are listed in the last column of Table 1). 

Considering the pair ^^^Cd-^^^Sn, the GT^^^ strength distribution calculated within 
the pn-dQRPA for the transition ^^^Sn— *^^^Sb is shown in Fig. la (a small imaginary part 
is added to the s-p potential). To compare the calculation results with the ^^^Sn(^He,t) 
experimental data of Ref. ^, five centroids of the energy, E^^i, and their strength Xi 
relative to the one of the GTR are evaluated (Fig. lb). The value g' = 0.77 allows to 
reproduce the experimental GTR energy in the calculation. The GT*^^-* strength distri- 
bution is almost insensitive to the gpp value {gpp = 1.0 is taken in the calculation). The 
GT*^"*"^ strength distribution for the transition ^^^Sn^^^^In is found more sensitive to gpp. 



0.47 corresponding to the lgg/2 



lgy2 transition 



Only one 1+ state with B^+\GT) 

into the ^^^In ground state, is found in the calculation within the interval E^ < 5 MeV. 
This weak transition is allowed due to the neutron pairing in ^^^Sn. In the ^^^Sn(d,^He) 
experiments four 1+ states in ^^^In were found within the interval E^ < 3 MeV with 
total strength Y.i Bi (GT) = 0.66 [3J. Population of the 1"*" states in ^^^In has also been 
studied in the ^i^Cd(p,n)-reaction ^2]. The result B^-\GT) = 0.26 ± 0.02 for excitation 
of the ^^^In ground state is only available now. Within the interval E^ < 3 MeV the 
calculated GT'-"-' strength distribution in ^^^In exhibits one 1^ state, corresponding to 
the back-spin-flip transition lg^/2 —>■ ^9q/2 i^^o the ^^®In ground state with the value 1.05 



/9/2 



B^ \GT). (for ^^®Sn-^^^Sb this transition is Pauli blocked). The 2z//5/3-decay amplitude 



for the decay ^^^Cd^^^^Sn can barely be evaluated within the pn-QRPA because the 
proton shell is closed in ^^^Sn. 

Coming to the pair ^^''Te-^^'^Xe, the value g' = 0.88 is found in the calculation by 
fitting the experimental GTR energy in ^^°I [9j. Then the amplitude M^^ (6) (or (9)) 
and its decomposition (10) are calculated, as a function oi gpp (Fig. 2). The corresponding 
experimental value {M^rpY^P = 0.03 MeV~^ [TU] can be reproduced in the calculation at 
dpp = 0.99. The 2z//?/3-decay strength function uj~^S^ \uj) is calculated for this value of 
gpp (Fig. 3a) along with the corresponding running sum Mo'rp^Lu) = J'^ uj'~ S^ \uj')duj' 
(Fig. 3b). Figs. 2 and 3 illustrate how the M^^ value for the decay ^3°Te->^^°Xe is 
formed. In particular, as one sees in Fig.3, the experimental studies of B^^\GT) are not 
always sufficient for understanding partial contributions to M^^. The reason is that the 
intermediate states having a relatively large excitation energy and very small B^^\GT) 
value (like the GTR) can nonetheless play essential role in formation of the 2z//?/3-decay 
amplitude. 

In conclusion, an isospin-selfconsistent version of the pn-cQRPA has been outlined and 
some its applications to description of charge-exchange excitations in open-shell spherical 
nuclei are presented. Although only general features of the low-energy strength distribu- 
tions can be described within the approach, it seems applicable to analysis of /3/5-decay 
observables. 
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Figure 1: The GT*^^) strength function for ^^^Sn-^^^Sb (a) and the relative (with respect 
to the GTR) strength of the low-energy l"*" peaks calculated within pn-cQRPA (b). The 
corresponding experimental data are taken from Ref. [S] 
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Figure 2: The calculated amplitude of ^'^°Te 2i//3/3-decay as a function of gpp. Decompo- 
sition of Eq. (10) is also shown. 
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Figure 3: The GT 2z//5/3-decay strength function (a) and the running sum (b) calculated 
for i30Te at gpp = 0.99 



